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We propose an Ansatz for universal conductance fluctua- 
tions in continuous dimensions from up to 4. The Ansatz 
agrees with known formulas for integer dimensions 1 , 2 and 
3, both for hard wall and periodic boundary conditions. 
The method is based solely on the knowledge of energy 
spectrum and standard assumptions. We also study nu- 
merically the conductance fluctuations in 4D Anderson 
model, depending on system size L and disorder W . We 
find a small plateau with a value diverging logarithmically 
with increasing L. Universality gets lost just in 4D. 

PACS numbers: 71.30.+h, 71.23.An, 72.15.Rn 

Disordered systems usually possess the metallic regime, 
where Ohm's law for mean conductance (g) works well at 
least for cubic samples and the distribution of g for various 
realizations of the same disorder is Gaussian with constant 
(disorder, mean free path l e and (g) independent) width, 
called universal conductance fluctuations (UCF). Higher 
cummulants of g should disappear as some power of l/(g) 
[1], though recent experiment on gold wire did not con- 
firm this in quasi-lD [2]. Recently [3] we analysed the 
statistical properties of conductance on bifractal lattices 
[4]. It became clear, that (g) and var g depend (besides 
the spectral dimension d s ) on lattice topology. Simply 
speaking, bifractals are no hypercubes. But we hope that 
by changing other parameters, say anisotropy, we can tune 
the systems to cubic-like. If the UCF of anisotropic bifrac- 
tals become those of non-integer dimensional hypercubes, 
other parameters of these systems may be comparable. 
The main goal of this work is to find a way to calculate 
the UCF also for non-integer dimensions. Rewriting (2) 
as a simple integral, we will propose its analytical contin- 
uation to real dimensions. 

Two (on first sight) different expressions were given for 
UCF in literature [5] and [6] . Let us comment this seeming 
ambiguity. In classical papers of Lee, Stone and Fukuyama 
[5] - Appendix, formulas forUCF in 3D (2D, ID) were 
given as a sum of three diagrammatic terms F a ,F and 
F c . They can be written in terms of convolutions, e. g. 
(the simple one- loop "Meeron" diagram): 



F a = 



L 4 



J n(r,r')n(r',r)drdr' 



(1) 



where II(r, r') is standard propagator in a box, F and F c 
contain 3 and 4 propagators (cyclically in one loop), re- 
spectively. In ID it is possible to show analytically, that 
F = —F a and F c — 3/4 F a , with some tricky cancelling of 
divergences, Ref. [7]. Numerical calculations show, that 
the same holds for 2D and 3D. Thus we arrive at the for- 
mula, given in Ref. [6] and representing 3/4 F a from Ref. 
[5]: 
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where the number of sums defines the dimensionality, hard 
wall boundary conditions are applied and we limit our- 
selves to cubes (L x — L y — L z ). This relation was de- 
rived using the technique of Hikami boxes. It gives in one 
dimension 

(9 2 )l D = ~ A C(4) = ^ (3) 

7T 4 15 

C(x) is Rieman's function. No such simple formulas were 
given for higher dimensions yet. It is also clear, that a four- 
fold sum, i. e. 4D case of Eq. (2) diverges and d = 4 plays 
a role of some kind " critical dimension" . According to the 
remark 6 of Ref. [8] , it is possible that also ergodicity gets 
lost at that point. 

Dealing with diffusive part of conductivity itself, the 
following formula for diagonal part of the propagator in 
3D was given [9] and [10]: 
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where Sg is the diffusive part of conductivity and a func- 
tion S(y) was specified as follows: 
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Here y — ir<fid(le/L) 2 , (fid is dimension specific, e. g. 1/3 
in 3D, but for our purposes insignificant constant. HW 
stands for hard wall and PBC for periodic boundary con- 
ditions in directions perpendicular to the current, flowing 
in the z direction. Later we will give an alternative mean- 
ing to the variable y. Contrary to Ref. [9] we include 
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zero modes in our Eq. (6), i x = and i y = 0, similarly to 
Ref. [11]. Our normalization of y differs slightly from that 
in Ref. [9]. Instead of complete elliptic integrals, we will 
make use of $3 from the family of Jacobi elliptic functions, 
defined as [12]: 



tf 3 («,g) = l + 2y q n cos(2nw) 
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where we set q — exp(— ny) and u = 0, making (7) related 
to modular elliptic functions. It is worth mention that all 
# functions fulfil the ID diffusion equation, 1?" + cd = 0. 

Let us recall the time t dependent retarded Green's func- 
tion (propagator): 

G(r, t, r', 0) = 0(i) exp(-iE n t)* n (r)K(/) (8) 

n 

where 9(t) is Heaviside step function, £"„ are eigenvalues 
and orthonormal complete eigenfunctions of the proper 
Hamiltonian and fi = 1. In our scaling, e. g. 3D HW box 
case, E n = n^ 2 . + i y + i 2 ). One can go over to spectral 
representation performing a Fourier transformation 
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G(t,t',E) = -i exp{iEt)G(r,t,r',0)dt 
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Inserting this into (1) and exploiting the orthonormality 
one gets the Meeron diagram value 



(e - E n y 
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as it was done in [5] and [6] at E = 0, see Eq. (2). The con- 
stant C, including appropriate combinatoric factor, will be 
specified later. 

We have an alternative way. Let us integrate out the 
space variables in advance 



Z(t) 



drdr'G(r,0,r',0)G(r,t,r',0) 



exp(-iE n t) 
(11) 

This simplifies the calculation, we will not have problems 
with the divergence of II(r, r) cx S(y) in more than one 
dimension, or with non-uniformly convergent, r-dependcnt 
series, Ref. [13]. Now we either first perform the Fourier 
transformation to get Z(E) = ^2(E — En)^ 1 and then 
apply a partial derivation, or we take the derivative of the 
Fourier integral to get 



Fa(0) 
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where y now means imaginary time y = it, thus going over 
from (zero potential within a box) Schrodinger equation to 
the diffusion equation, or from Fourier to Laplace transfor- 
mation. The quantity Z(y) strongly resembles a partition 
function, with y playing the role of inverse temperature. 



Imry [14] also stated that variance of g over random 
matrix ensembles can be calculated by a similar integral 
formula. Our version of (1) and thus (2) with proper con- 
stant now reads: 
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Noticing that Z(y) = J2cxp(—E n y) — —S'(y)/ir, after 
per parts we get an interesting, though for calculations 
not quite practical formula 
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S(y)dy. 
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Specifying Z(y) in (13) with help of $3 function (7) we 
can rewrite (2) for small integer dimensions d as follows: 
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y[M0,q) - 1] [#3(0,<z) + If-'dy. 



(15) 

This is easy to verify: each term with $3 creates one sum- 
mation index in (5) for dimensions d = 1, 2 and 3. After 
per parts and integrating term by term we get (2). 
In ID case, one can make use of the formula [12] 
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If we put y — x 2 in (15), specifying s — 4 yields (3). 
Similarly for PBC 
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(9 2 )c BC = ^J o y[U0,l)-mU0,q 4 )} d - 1 dy. (17) 

Note that q 4 = exp(— Any). Now we can calculate the UCF 
with high numerical precision for small integer dimensions, 
see the Table 1. But the main generalization consists in 
regarding (15) and (17) as an analytical continuation in d, 
i. e. being valid also for non-integer dimensions [3] and 
even for those ones below ID, where it is hard to imagine 
any "perpendicular" direction. 

The Equations (15) and (17) represent the completed 
Ansatz. It reproduces the UCF for any integer dimension, 
i. c. any number of sums in (2), which should make the 
continuation unambiguous. But this argument is rather 
naive, as the sums diverge for d > 4; some cut-off (e. 
g. as in [9]) would be necessary. Note that the integrals 
(15) and (17) also diverge for d > A. The dependence 
of UCF on dimension is plotted in Fig. 1. The common 
statement, that this dependence is weak [6], is true in the 
region < d < 3. 

We shall now comment the limiting cases. For d = 
one gets sound values (see Table 1) both slightly below 
the universal constant 1/8, which corresponds to ballistic 
transport in quantum dots [15]. For d approaching 4, d = 
A — e the leading contribution to integrals (15) and (17) 
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comes from small y values and we can make use of the 
known alternative expansion: 

1 00 ( ^2 

Mu,eM~™ 2 )) = l E ^V J ^? L ] (18) 
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For small x — y/y and u = only the n = term is 
important (see Ref. [9]) 

d 3 (0,exp(-7rx 2 ))^- (19) 

x 

and we get 

<S 2 >c = ~ + 0(l) (20) 

both for HW and PBC. These expansions remain valid 
even for anti-periodic BC, if one simply replaces the sec- 
ond $3 function in (17) and that one in (19) by i?2- The 
independence of the leading divergent term on boundary 
conditions seems to support our Ansatz. Similar state- 
ment was made in 3D for small y expansion of S3 (y) , whose 
leading divergence also proved BC independent, see [9]. 

Unfortunately we cannot use formulas (15) and (17) to 
get the exact i-dependence of 4D conductance fluctua- 
tions directly. Introducing 1/L cut-off, e. g. by changing 
the lower limit of integration to small non-zero value, cor- 
responds to making upper limits for summation indices in 
(2) finite. This would influence the propagator used in 
one-loop diagrams so that the delta function on rhs of dif- 
fusion equation in spectral representation (Refs. [5], [6]) 
would become just a peak of finite, L dependent width. 
Thus the summation of diagrams, mentioned in the begin- 
ning, would be no more simple. Anyway, the too simple 
cut-off would give logarithmic L dependence in 4D. 

In [16] numerical calculations of UCF for 2D and 3D An- 
derson model were compared to theoretical values, com- 
puted also numerically with cca 3 digit precision from more 
complicated version of (2), given in [5] and its PBC coun- 
terpart. A plateau of var g = (g 2 ) — (g) 2 = (g 2 ) c as a 
function of disorder W was found with values reasonably 
close to theoretical predictions. The disorder had to be 
large enough to overcome the ballistic peak region, but 
not too close to the metal - insulator transition (if > 2). 
The plateau grew broader with increasing system size L. 

We performed these calculations for 4D Anderson model 
with HW BC. Up to now only the region of metal - insula- 
tor transition was addressed in 4D, Refs. [3], [4] and [17]. 
In diffusive regime, Ohm's law for cubic samples applies 
in the sense of leading contribution in L: 

(g) = oL d - 2 (21) 

This serves as a test of diffusive kind transport, we expect 
(g) <x L 2 , see Fig. 2. Here we find, that this regime is well 
pronounced around W = 11. 

Results for var g are shown in Fig. 3. A small plateau 
appears for larger L and disorders cca within 10 < W < 
11. Typical statistical ensembles are N 8tat = 10 5 for L < 



6, w 20000 for L = 7 and w 3000 for L = 8. The value of 
the plateau diverges logarithmically (Fig. 4). 

We conclude that a simple conjecture enables a high pre- 
cision calculation of UCF in any real dimension < d < 4, 
with reasonable behavior in both limiting regions. Numer- 
ical calculations of 4D Anderson model show logarithmic 
divergence of CF with increasing system size. 
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Figure 1: Universal conductance fluctuations as a func- 



tion of dimension for HW (open symbols) and PBC (full 
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L Figure 4: Conductance fluctuations in 4D as a function of 

L. The slope varies from 0.36 to 0.39. 

Figure 2: Mean conductance in 4D as a function of system 
size. Fitting parameter n from (g) oc L n (for L > 5) shown 
in the Figure. 
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